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Abstract 

We examine dijet production at large rapidity intervals at Tevatron energies by compar- 
ing an exact O(o%) calculation with the BFKL approximation, which resums the leading 
powers of the rapidity interval y to all orders in a s . We analyze the dependence of the 
exact 0(a^) calculation on the jet cone-size as a function of y, and use this cross section 
to define an "effective rapidity" y which reduces the error that the large-?/ approxima- 
tion induces on the kinematics. Using y in the BFKL resummation, we reexamine jet 
production at large transverse momenta and the transverse momentum decorrelation of 
the tagging jets. We find less dramatic, but still significant, effects than found previously 
using the large-?/ approximation. 
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1 Introduction 

The state-of-the-art in jet physics at hadron colliders is described by next-to-leading- 
order QCD parton-level calculations. These consist of the 0(a 2 s ) + 0(oP s ) one-loop 
2 — > 2 parton scattering, combined together with the 0(a^) tree-level 2 — > 3 parton 
scattering]]]]. These perturbative calculations describe the hard part of the scattering, 
while the nonperturbative effects are factorized into the parton structure functions. An 
advantage of going to next-to-leading order is that it reduces the dependence on the 
arbitrary scale associated with this factorization. In addition, the inclusion of a third 
final-state parton allows a more detailed description of the jet structure. These next-to- 
leading order calculations @, ||, || appear to be in very good agreement with the one- 
and two-jet inclusive distributions obtained from the data of the CDF experiment at the 
Fermilab Tevatron Collider ||, |J. 

Despite these successes, it is possible to imagine kinematic configurations where this 
fixed-order analysis is inadequate, even though the underlying process is still perturba- 
tive in the usual sense. This could occur when the cross section contains logarithms of 
large ratios of kinematic invariants. Typical invariants are the hadron-hadron center- 
of-mass energy y/s, the parton-parton center-of-mass energy Vs = ^xaXbs, where xa 
and xb are the momentum fractions of the partons originating the hard scattering, and 
the momentum transfer Q, which is of the order of the transverse momentum of the jets 
produced in the hard scattering. Large logarithms will appear when yfs ^> Q, in the 
semihard region of kinematic phase space. We can then avoid complications associated 
with the small-x behavior of the structure functions by requiring that the parton mo- 



mentum fractions, xa and xb, are sufficiently large, as originally suggested by Mueller 
and NaveletQ. In this case the large logarithms, ln(s/Q 2 ), factorize entirely into the 
partonic subprocess cross section. These logarithms, which are of the size of the rapidity 
interval in the scattering process, can be resummed by using the techniques of Balitsky, 
Fadin, Kuraev, and Lipatov (BFKL)||. For jet events with a large rapidity interval, the 
amplitudes are dominated by contributions from multiple gluons that uniformly fill the 
interval between the two extreme jets. The BFKL theory systematically resums these 
leading powers in the rapidity interval, including both real and virtual gluon corrections. 

In a previous paper[[|] we showed how to analyze dijet production experimentally so 
that it most closely resembles the configuration assumed in the BFKL theory. The main 
difference from the standard hadronic jet analysis is that the jets are ordered first by 
their rapidity rather than by their energy. Thus, we look at all the jets in the event that 
are above a transverse momentum cutoff p± m i n , using some jet-definition algorithm, and 
rank them by their rapidity. We then tag the two jets with the largest and smallest 
rapidity and observe the distributions as a function of these two tagging jets. The cross 
section is inclusive so that the distributions are affected by the hadronic activity in 
the rapidity interval y between the tagging jets, whether or not these hadrons pass the 
jet-selection criteria. We will refer to these hadrons in the rapidity interval as minijets. 

In ref. we showed that the exponential enhancement with the rapidity interval y in 
dijet production at fixed xa and xb , which was originally suggested as a signature of the 
BFKL minijets by Mueller and Navelet, is highly suppressed by the parton distribution 
functions at Tevatron energies. However, other observables such as the jet transverse 
momentum distribution, and the jet-jet correlations in p± and azimuthal angle are signif- 



icantly affected by the minijets. For instance, the transverse momentum distribution was 
considerably enhanced at large p± and large y This enhancement has been seen in the 
CDF data||, although the data analysis differs somewhat from that needed to compare 
directly with the BFKL resummation. In addition, there should be some dependence on 
the cutoff p±_mi n , and we shall see here that the p± distributions should differ for the two 
tagging jets if they are not symmetrically placed around zero rapidity. Finally, we saw 
that the correlation in transverse momentum and azimuthal angle of the tagging jets 
is not a leading feature of the expansion in the rapidity interval. Accordingly, it fades 
away as the rapidity interval increases. The decorrelation in azimuthal angle has also 
been noted by Stir ling fT0|j. 

For the most part these effects at large y can be easily understood in terms of the 
sharing of the total p± by the additional minijets, which the BFKL resummation auto- 
matically includes. In a fixed-order calculation the effects arise first at 0(a^), where the 
introduction of a third final-state parton removes some of the correlations inherent at 
lowest order. Thus, it is interesting to compare the BFKL resummation directly with the 
0(a^) result. This was done for the case of the phi distribution in ref. [|Kj. Here we shall 
extend this comparison, which can be approached in two different ways. First, one can 
truncate the BFKL solution at 0(a^) and compare this with the exact 0{a^ s ) result. The 
truncated 0(a^) BFKL cross section is just the large y limit of the complete 0(a^) cross 
section, and so this comparison will indicate how good the leading log approximation 
is at the different experiments. Second, one can compare the full BFKL solution with 
its truncation to 0(oP s ). This should give an indication of the size of the contributions 
from higher orders in a s . In addition it will help in isolating the distributions which are 



sensitive to the resummation effects beyond next-to-leading order. 

The remainder of this paper is as follows. In section 2 we discuss the BFKL solution 
and its truncation to 0(oP s ). For illustrative purposes we also show how this second 
cross section is obtained from the 0(a^) tree-level 2 — ► 3 cross section in the limit of 
large y. In section 3 we discuss the infrared singularities that occur at 0(oP s ) in order 
to further elucidate the approximations used in the BFKL analysis. We also observe the 
dependence on the jet cone-size as a function of y. In section 4 we compare the BFKL 
solution truncated to 0(af) with the exact 2 — > 3 0(a^) cross section. We isolate the 
primary contribution to the discrepancy between the two as arising from approximations 
used in the parton distribution functions. In section 5 we introduce an "effective rapidity" 
variable to use in the BFKL formula which accounts for this difference. In particular, 
the effective rapidity is defined so that the truncation of the BFKL solution to 0(a^) 
equals the exact 2 — > 3 0(a^) cross section. We use this modified BFKL cross section 
to study the various dijet distributions, which show less dramatic, but still noticeable, 
effects than in ref. Q. In section 6 we present our conclusions. 

2 The minijet resummation and 0(a^) cross section 
at large rapidities 

We are interested in the semi-inclusive production of two jets in hard QCD scattering. 
For definiteness we will consider the scattering process PaPb —> 3\3i + X such as at the 
Tevatron, but the same analysis can also be applied to photoproduction at HERA[|i~l"|]. 
We describe the two partonic tagging jets by their transverse momenta and rapidities 



(pi±,yi) and (j>2±,y2), where we always take y\ > y<i- For large rapidity intervals, 
y = yi — 2/2, the cross section for this process can be written 



^2 _ „0„0 f (Ji 2\ f / ,.2\ d&gg 



x° A x° B f eS (x° A , v 2 )f eS (x% f) "? 2 3 j- , (1) 



where the parton momentum fractions are dominated by the contribution from the two 
tagging jets 

o Pi±e Vl 

Xa ~ ~^r 

A ~ ^ (2) 

and n is the factorization/renormalization scale. In this limit the amplitude is dominated 
by gg, qg, and qq scattering diagrams with gluon-exchange in the t-channel. The relative 
magnitude of the different subprocesses is fixed by the color strength of the respective 
jet-production vertices, so it suffices to consider only gg scattering and to include the 



other subprocesses by means of the effective parton distribution function fT2| 



Q 

f cS (x, n 2 ) = G(x, n 2 ) + — - £}[Q/(x, n 2 ) + Q f (x, /i 2 )] , (3) 

Ca f 

In (|3p the sum is over the quark flavors, C A = N c = 3 is the Casimir operator of the 
adjoint representation and Cf = (N% — 1)/2N C = 4/3 is the one of the fundamental 
representation. 

The higher-order corrections to the gg subprocess cross section in ([!]) can be expressed 
via the solution of the BFKL equation^], which is an all-order resummation in a s of the 
leading powers of the rapidity interval 

d(T 99 _ C A a s S^ JnU-TT) f°° j,.„u(n,i/)y _ „ / ,. i„ PlX 



Y, e in{ *- w) f" dve"^ y cos ( v In ^ ) , (4) 



dp 2 ± dp 2 2± d(j) ^Pi ± pl ± n J ® V p\l 



with 



u[n, v) 



H^-^(i)_Re^H+l + il/ ) 



(5) 



n 2 

and i[) the logarithmic derivative of the Gamma function. Eq. (Q) can be expanded order 
by order in a s and compared with a fixed-order calculation of dijet production at the 
same order of a s , in the large rapidity limit |T0|. By expanding the exponential in Q) to 
zeroeth order in a s we obtain the tree-level large-?/ cross section 
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At 0(a s ) in the exponential we obtain 
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When pu_ ^ p 2 ± we can integrate out v and sum explicitly over n, by using the integral 
representation of the ip function 



V(z) 



ax 

o 1 — X 



with the Euler constant 7 = .577215 . . .. We obtain 

dai 1 ) C\ot 2 
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C A ot s y 



(8) 



dpf ± dpi ± dcj) 4:7Tpf ± p?. x pj ± + p( ± + 2p 1± p 2± cos . 



(9) 



When pu_ ~ p 2 ±_ the cross section (|7]) is dominated by configurations where the third 
parton is soft. These infrared singularities are regulated by the BFKL solution, which 
includes both real and virtual corrections. To see how this occurs we can integrate (0) 



over £>2_i_ m the interval defined by \pfj_ — p|j_ | < epf ± , for e sufficiently small. The 
integrals over v and x can be performed, and the series in n can be summed, giving a 
finite answer: 

S(</>-7r)—^ylne, (10) 



dag) «CW 2C A a s 

dp\ ± d<j> 2 P \ ± 0{9 7T) TX 



up to terms of 0(e). For sufficiently small e this configuration is indistinguishable from 
a configuration with only two final state partons. Thus, for \pfj_ — p\^ \ < epfj_ we can 
write the cross section to 0(a 3 ) as 

{ s 5(p 2 lx -pl x )5(<f)-Tr) [1 + ylne . (11) 



dp\ ± dp2 ± d(f) 2pf± v ' ri " L l ' ZJ ~ J ' V n 

For |pjj_ — p2_|_| > epfj_ we can use equation (H). Combining these two formulae (|ll|) 
and (0), the dependence on the unphysical variable e will vanish in any inclusive process 
integrated over a range of momenta, for small enough e. 

It is also informative to derive the 2 — > 3 cross section (|S|) by taking the large-rapidity 
limit of the 0(a s ) real corrections to gluon-gluon scattering, computed in the conven- 
tional way||13||. First, note that in the large rapidity limit the leading contribution is 
given by the gg — > ggg subprocess, the diagrams with final-state quarks being sublead- 
ing. The squared scattering amplitude for this, summed (averaged) over final (initial) 
colors and helicities, is given in ref. |14| as 



\Mf = 4(na s C A rY.4 E Sa s ' s Sa , (12) 

i>3 [4,1,2,3,*]' S A\Sl2S2ZS ZB SAB 

with i,j = A, 1,2, 3, B, and with the second sum over the noncyclic permutations of the 

set [A,1,2,3,B]. The kinematic invariants are defined here by 

it 
sab = s = x A x B s = ^ Pi,±Pj,±e Vt ~ Vj 



s A ,i = -J2PiM± e ~ ivi ~ K) ( 13 ) 

n 
3=1 



s 



t.l 



2Pi,±Pj,± [cosh(j/j - Vj) - cos(0i - 0j)] . 



In the large rapidity limit eq. ([E|) becomes | 



,s 2 



\M\ 2 = 128 (na s C A ) 3 2 \ 2 . (14) 

The large rapidity limit of the phase space for three-particle production is 

^--h d -B(uMh r <P-^ (i5) 

where we have fixed the rapidities of the two gluons at largest and smallest rapidity and y^ 
is the rapidity of the third gluon. Using the amplitude (|14D , the phase space (|i~5[) and the 
appropriate flux factor, we obtain for the dijet production cross section da g g/dpl ± dp2i_d<t> 
the same expression as the one obtained in eq. (Q) via the expansion of the BFKL solution. 
Note that the overall factor of y in eq. @ comes from the integration of the rapidity of 
the third final-state gluon over the interval spanned by the tagging gluons. 



3 Collinear singularities in dijet production 

Before entering the details of a comparison between the BFKL resummation and the 
complete 0(af ) result, it is useful to consider the infrared singularities that occur in the 
two approximations to dijet production. These singularities depend on the geometry of 
the event and thus are sensitive to how dijet production is defined. In the next-to-leading 

8 



order (NLO) approximation the phase space integration over the kinematic variables of 
the third (unresolved) parton generates infrared singularities when it becomes soft and 
collinear singularities when it becomes collinear either with the initial-state partons or 
with the other final-state partons. In contrast, the BFKL approximation only contains 
the effects of soft gluons. Because of the explicit rapidity-ordering, the unresolved partons 
are bound to lie within the rapidity interval defined by the tagging jets, and so can never 
become collinear with the initial-state partons^. In addition, in the large rapidity limit 
the region of phase space where the third parton becomes collinear with either of the 
other final-state partons gives a nonleading contribution. 

Let us examine this last point further. First we recall the standard definition of a jet as 
consisting of all the partons lying within a circle of radius R cut = [(yj — y) 2 + {4>j — <fi) 2 ] l ' 2 
from the jet axis (yj, <pj) in the plane defined by rapidity and azimuthal angle (i.e., the 
Lego plot). Then we look at the dependence on the jet cone size R cut of the 0(a^) 
corrections to dijet production. Namely, we consider in Fig. 1 the contribution of the 
three-jet configurations, computed through the 2—^3 parton amplitudes, to the cross 
section 

da fvi , f , , ^-^ „ 0s „ 0s d&ij 



/ dy 3 d<f)Y, x AX B fi/A{xA,^ 2 )fj/B{xB,^ 2 

Jw> J 



dydydp 1 _ L dp 2 _ L hi J ^ dp 1± dp 2 ±dy 3 d(p 

(16) 

as a function of the jet cone size R cu t at different values of the rapidity interval y. 

fi(J) = QiQ->G labels the distribution function of the parton species and flavor i(j) = 

2 Of course, it is always possible for a parton or jet to be produced outside of the rapidity interval of 
the tagging jets, but to go undetected because it has too small pi± or too large \yt\. In the context of 
the BFKL approximation this jet would be considered part of the evolution of the structure functions, 
and would not be included in the calculation of the parton subprocess. 



q, q, g inside hadron A(B). We include all parton subpro cesses [fX3|] , and use the exact 

values of the parton momentum fractions 

Pi±e m +p 2 ±e m + p 3 ±e V3 



x A 



V~s 



p 1± e yi +p 2 ±e y2 +p 3 ±e vz 
x B = 1= • (17) 

As in our previous analysis we set the rapidity boost y = (yi + y 2 )/2 = 0, since we are 
mainly interested in the behavior of the parton subprocess, which does not depend on y. 
We fix the transverse momenta of the tagging jets at different values, p±±= 20 GeV and 
P2±= 50 GeV, in order to ensure that the third parton cannot become soft. Furthermore, 
collinear configurations where the distance R between two of the partons on the Lego plot 
is smaller than R cut are discarded, because the resulting "combined" jet would be back- 
to-back with the remaining parton and have the same p± as it. Finally, in the spirit of 
the BFKL approximation, we restrict the rapidity of the third parton to lie between the 
two tagging jets. Note that this last approximation differs from the full NLO treatment 
which includes initial-state collinear radiation and uses the Altarelli-Parisi subtraction to 
handle the divergences. Therefore, in this calculation and in all BFKL calculations, the 
parton density functions should be treated at leading order (LO) and will have all the 
associated LO factorization-scale dependence. In all of our plots we use the LO CTEQ 
parton distribution functions JTJJ] with the renormalization and factorization scales set to 
the geometric mean of the transverse momenta of the tagging jets /z 2 = pi±p 2 j_. As Fig. 
1 shows, the dependence of the dijet production on the jet cone size R cut decreases as y 
is increased. At y = 6 the tagging jets are so widely separated in rapidity that there is 
almost no dependence on R cut at all. This confirms that the phase space region where 
two final-state partons are collinear is subleading at large rapidities, as suggested by the 
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BFKL theory. 

4 The BFKL and the 0(af) corrections to the p± 
distributions 

In section 2 we have shown that the truncated 0(a^) BFKL contribution to dijet pro- 
duction is the large rapidity limit of the complete 0(a^) corrections. Here we want to 
see how well this works in practice in the jet-jet p± correlation case examined in ref.§|. 
In Fig. 2 we consider the contribution of the three-jet amplitudes to the p± correlation 
da/dydydpi^dp2±, plotted as a function of the transverse momentum pi±, at a fixed 
value of P2±= 50 GeV and at y = 2 and 6. The customary value for the jet cone size 
Rcut= 0.7 has been used. As in fig. 1, configurations where the distance R between two 
partons on the Lego plot is smaller than R cut are discarded, since they would be counted 
as a two-jet event with pi± = p 2 ±. We go from the exact configurations to the large-?/ 
approximation to them in three steps. The dashed curves are computed through the 
exact 2— >3 parton amplitudes |13| and kinematics ([17]); the dotted curves through the 



large-?/ parton amplitudes (|TJ]) and the exact kinematics (0); the solid curves through 
the large-?/ parton amplitudes (|T4]) and kinematics (0). As the plots show, the error in 
using the large-?/ approximation grows with the imbalance in transverse momentum of 
the tagging jets. While at small ?/'s the error is distributed between the approximation 
on the amplitudes and the one on the parton distribution functions, at large ?/'s most of 
the error comes in using the large-?/ approximation in the parton distribution functions. 
We have also made this comparison for the larger rapidity intervals obtainable at the 



11 



CERN Large Hadron Collider (LHC) and have found that this discrepancy, although 
smaller, is not insignificant. 

This discrepancy can be understood by recalling that the rapidity 2/3 of the third jet is 
integrated over the full range of the interval from y 2 to y\. If we neglect its contribution 
to the momentum fractions as in (|2|), this just multiplies the cross section by a factor of 
y as in eq. (^). However, eq. (0) can be a bad approximation to the exact kinematics 
(|17|) over much of the integral if p%± is not small. For y% near the extremes, using the 
exact kinematics in the parton distribution functions produces a large suppression, so 
that the "effective" rapidity range of 2/3 is reduced substantially. The truncated 0(a^) 
BFKL (solid) curve neglects this effect, and so it overestimates the cross section. Note, 
however, that near p\± = p 2 ± the transverse momentum of the third parton is small, so 
its contribution to the rr's in (O) can be safely neglected. 



5 An effective rapidity interval 

We have just seen in the previous section that the large-?/ approximation used in the 
BFKL resummation seriously overestimates the cross section when the two tagging jets 
have unequal energies, even for rapidity intervals as large as y = 6. This occurs because 
the large-?/ cross section (|j|) assumes that the third (minijet) parton can be produced 
anywhere within the rapidity interval [1/2,1/1] with equal probability, whereas in the full 
2^3 cross section the probability is highly suppressed by the structure functions when 
the third jet strays too far from the center of this interval. In this section we will attempt 
to fix this problem by including this suppression effect directly into the BFKL equation. 
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The BFKL solution as presented in equation (H) is an all orders resummation in 
a s y where y is the kinematic rapidity interval defined by the tagging jets. However, 
to be completely precise, the rapidity variable which is resummed by BFKL is only 
defined up to transformations y — > y + X where X is subleading at large rapidities. 
This is analogous to a change in the scale Q 2 in the standard a s ln(Q 2 ) resummation, 
and in the same way that one chooses a physical value of Q 2 to lessen the effects of 
nonleading terms in this resummation, one can choose a more physical rapidity variable 
in the BFKL resummation. From the results of the previous section we see that a better 
rapidity variable would be one that reflects the range in rapidity spanned by the minijets, 
which is typically less than the kinematic rapidity interval y. Let us define an "effective 
rapidity" y(n,p 1± ,p 2 ±,y,y) by 



dcj) cos(n<f>) (da / dydydpi±dp 2 ±d(f)) 

y = y - r , (is) 

dcf) cos(n<j)) (da /dydydpi±dp2±d(p) 



where n is the fourier series index of eq. (0). The cross section in the numerator is that 



of eq. (|16| ) and is computed using the exact kinematics (|17|), while the cross section in 
the denominator is that of eq. ([I]) and is computed using the large-?/ kinematics (|2]) . The 
denominator can easily be computed analytically using the large-?/ solution (||). Note 
that y is defined so that if we replace y — > y in the BFKL solution (|j) and truncate to 
0(a z s ) we recover the exact 2^3 cross section. Also note that asymptotically for large 
y, the difference y—y is nonleading. Thus, we can use y in (|) and obtain a quantitatively 
more reliable solution. 

In Fig. 3 we plot y as a function of pi± for n = 0, p 2 i_ = 50 GeV, y = 0, and 
y = 2,3, 4, 5, 6 with R cut = 0.7. For y = 2 and 6 this is just given by the ratio of the 
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dashed curves to the solid curves in Fig. 2, multiplied by y. Near pu_ = p 2 j_, y approaches 
the kinematic rapidity y, especially for large y, but it falls quickly as the two jets move 
apart in transverse momentum. Thus, we would expect the effects of the resummation 
to be most important near pi± = p 2 ±. We can see this clearly in the jet-jet p± correlation 
plot of Fig. 4, where we compare the exact 2 — > 3 0(a^) cross section (dashes) with the 
BFKL resummation using y (solid). As in Fig. 2 we fix p 2 ± = 50 GeV, y = 0, and plot 
as a function of pi± for y = 2 and 6, using a jet cone size of R cut = 0.7. In this and 
in the plots that follow, y is fixed at n = by the integration over the azimuthal angle 
0. As expected, for y = 2 there is little difference between the two approximations. For 
y = 6 the higher orders of the BFKL resummation are very important near pi± = p 2 ±, 
but have less of an effect when the jets are mismatched in transverse momentum. 

Next, in Fig. 5 we plot the transverse momentum distribution da j dy dy dpi± of jet 1 
at y = 4 and 6 and y = 0. The dashed curves are computed through the exact Born- level 
matrix elements and the solid curves using y in the BFKL resummation and two different 
cutoffs for jet 2 transverse momentum, p 2 ± m i n = 10 GeV and 20 GeV. The effect of the 
radiative corrections is to release the p± distribution from the Born-level requirement 
that pu_ = p 2 j_. At a fixed order of a s and for a given pu_, jet 2 and the minijets 
try to have p± as small as possible in order to minimize their contribution to the x's 



(17) and thus maximize the value of the parton distribution functions. They are only 
constrained by the overall transverse momentum conservation. Eventually when the 
number of minijets is very high, or virtually infinite as in the BFKL resummation, the 
smallest value of p 2 ± is not fixed any more by transverse momentum conservation but 
by the minimum p± experimentally attainable ||. This explains the strong sensitivity of 
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the curves in Fig. 5 to the value of P2±mm0- The crossing of the curves with j) 2 imm = 20 
GeV below the Born curves for small values of p±± is a further manifestation of this 
decorrelation in transverse momentum. For pi± ~ 20 GeV, jet 2 very often radiates 
away enough energy so that p2_i_ < P2±min- Thus, events that would have been kept at 
the Born level are now discarded when higher orders are included. Finally, for sake of 
comparison, we replot in the dotted curves the p± distribution computed at y = 4 using 
y in the BFKL resummation, with p 2 imm = 10 GeV and 20 GeV |J. We note that, using 
the effective rapidity y rather than the kinematic rapidity y in the BFKL resummation, 
the enhancement in the p± distribution is considerably reduced. 

To examine further the kinematic effects on dijet production we plot in Fig. 6 the 
p± distribution at \y\ = 2 and y = 4 as a function of the transverse momentum of 
the jets at largest (\yj\ = 4) and smallest (\y c \ = 0) absolute rapidities, which we call 
the forward jet and the central jet, respectively. As in Fig. 5 we use two different 
cutoffs p_L m in = 10 GeV and 20 GeV for the transverse momentum integrated out. 
Since we only have changed y with respect to Fig. 5 the contribution of the parton 
subprocess to the p± distribution will be the same as in Fig. 5, but the contribution of 
the parton distribution functions will change. The dashed curve is computed through 
the exact Born-level matrix elements and the dotted and solid curves using y in the 
BFKL resummation. At the Born level the transverse momentum distributions of the 
forward and central jets are the same since pj± = p c ±, but they become starkly different 
when higher orders are included. The dotted curves are the distributions of the forward 



3 On the other hand, light-cone momentum conservation requires that pi± < Y / sexp(— y\), so that 
at large rapidity intervals the maximum value of p\± does not appreciably change going from the Born 
level to higher order corrections. 
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jet. They do not appreciably differ from the Born-level curve because the upper bound 
Pf± < V^ ex P{~\yf\) from light-cone momentum conservation is very restrictive and 
the phase space of the forward jet does not basically change going from LO to higher 
order corrections. On the other hand the solid curves, which represent the distributions 
of the central jet, show a huge enhancement and very different slopes as compared to 
the Born level. To understand this, we have to look at the kinematics of the higher 
order corrections. At 0{a z s ) it is possible to let p c j_ grow by taking the third jet at 
small absolute rapidity and almost back-to-back with the central jet, p^\_ ~ — p c ± and 
y 3 ~ 0, while taking the forward jet as slim as possible, pf± ~ P± m in- Then from light- 
cone momentum conservation p c ± max ~ (y/s — p± m in ex P(|y/|)) /2, which is much bigger 
than the Born-level upper bound. This picture generalizes to higher orders with the 
transverse momentum of the minijets balancing that of the central jet, while the forward 
jet is produced near p± m i n . In order to test this picture we consider in the dot-dashed 
curve the transverse momentum distribution for the Born-level production of two jets in 
the central region \y\ = 0.5 and y—1. The slope of this curve is similar to that of the 
solid curves, suggesting that the p c ± distribution is dominated by events with two central 
jets and one soft forward jet, where we tag on the forward jet and one of the central 
jets. In the BFKL resummation the second hard central jet is presumably replaced by 
a succession of minijets. Thus, it is misleading to compare the BFKL resummation to 
the Born level (dashed) for the p c ± distribution, since it is dominated by events with 
more than two jets. A better comparison in this case would be between the BFKL 
resummation and a NLO calculation. 
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6 Conclusions 

In this paper we have attempted to better grasp the range of validity of the BFKL 
resummation by comparing it with an exact 2—^3 0(a 3 s ) calculation. We saw that the 
dependence on the jet cone-size quickly becomes insignificant for moderate rapidities, just 
as is assumed in the BFKL formalism. However, the approximation to the kinematics 
in the parton momentum fractions causes a serious error in the BFKL predictions when 
the tagging jets are not back-to-back in p± and 0. In order to account for this error we 
introduced an effective rapidity y which restricts the phase space of the minijets in such 
a way that the truncation of the BFKL resummation to 0(a 3 s ) agrees with the exact 
2^3 0(a 3 ) calculation. 

Using the BFKL resummation with the effective rapidity y we have seen that the 
effects on the p± spectrum are not as dramatic as we had predicted previously. The diffi- 
culty in detecting these deviations from the Born-level computation are compounded by 
renormalization/factorization scale ambiguities, which are at least as problematic here 
as at the Born level. Because of the two scales defined by the tagging jets we could 
let /i 2 = Cpi±p2±, Cmax(p 2 ± ,p2_i_)) or some other choice where C is some constant of 
order one. Because of the relatively small deviations of the BFKL resummation with 
the effective rapidity y from the Born-level calculation, and the sizeable renormaliza- 
tion/factorization scale ambiguities in the BFKL approximation, we expect that a com- 
plete NLO calculation could give a more reliable estimate to the rawp^ spectra. However, 
we do note that much of the uncertainties due to the renormalization/factorization scale 
drop out in the ratios of cross sections so that, for instance, the predictions of the ratio of 
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the pij_ spectra with different P2±min cutoffs are probably reliable to ~ f 5%. In addition, 
the large corrections in the p± correlation and <fi correlation plots suggest that an NLO 
calculation would be inadequate in these cases, and that the BFKL resummation using 
the effective rapidity y should do a better job here. 
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Fig. 1: Dijet production as a function of the jet cone size R cu t, at pi r 
P2,±_ = 50 GeV, y = and y = 2, 3, 4, 5 and 6. 
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Fig. 2: p± distribution of jet 1 with the jet 2 transverse momentum fixed at 50 GeV, 
at y = and at y = 2 and 6. The jet cone size is fixed at 0.7. The dashed curves are 
computed through the exact 2— »3 parton amplitudes and kinematics; the dotted curves 
through the large-?/ parton amplitudes and the exact kinematics; the solid curves through 
the large-!/ parton amplitudes and kinematics. 
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Fig. 3: y as a function of p±± for fixed p2± = 50 GeV and n — 0, y — 0. The curves from 
bottom to top are for y — 2, 3, 4, 5, and 6. The jet cone size is fixed at 0.7. The upper 
end of the curve y = 6 is cut off because pi± reaches the kinematic limit. 
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Fig 4: p_|_ distribution of jet 1 with the jet 2 transverse momentum fixed at 50 GeV, 
at y = and at y = 2 and 6. The jet cone size is fixed at 0.7. The dashed curves are 
computed through the exact 2— >3 parton amplitudes and kinematics; the solid curves 
are computed from the full BFKL solution using y. 
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Fig. 5: p± distribution of jet 1 at y — and y = 4 and 6. The dashed curves are the 
exact Born-level p± distributions. The solid curves are the p± distributions computed 
using y in the BFKL resummation, with two different cutoffs for jet 2, the upper curve 
with P2±min— 10 GeV and the lower curve with P2± m in— 20 GeV. The dotted curves are 
the p± distributions at y = 4 only, computed using the kinematic rapidity y in the BFKL 
resummation, with p 2 imin= 10 GeV and 20 GeV. 
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Fig. 6: p± distribution at \y\ = 2 and y — 4. The dashed curve is the exact Born- 
level p± distribution. The dotted (solid) curves are the p± distributions of the forward 
(central) jet computed using y in the BFKL resummation, with two different cutoffs for 
the transverse momentum integrated out, the upper curve with p± m in = 10 GeV and 
the lower curve with p± m i n = 20 GeV . The dot-dashed curve is the exact Born-level p± 
distribution at \y\ =0.5 and y — 1. 
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